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This  study  was  conducted  to  investigate  radiation  and  chemical  reaction 
effects  on  unsteady  MHD  boundary  layer  flow  in  porous  media  over  an 
exponentially  stretching  surface.  In  the  study,  boundary  layer  flow,  heat 
and  mass  transfer  of  incompressible  viscous  unsteady  MHD  fluid  over 
porous  stretching  vertical  surface  in  the  presence  of  thermal  radiation 
was  investigated.  A  similarity  transformation  was  used  to  reduce  the 
governing  system  of  partial  differential  equations  (PDEs)  to  a  set  of 
nonlinear  Ordinary  Differential  Equations  (ODEs)  which  were  solved 
numerically  using  computer  embedded  fourth-order  Runge-Kutta  method 
with  shooting  technique.  The  results  showed  that  an  increase  in  radiation 
produced  a  rise  on  the  velocity,  temperature  and  concentration  profiles. 
Skin  friction,  Nusselt  and  Sherwood  numbers  increased  with  increasing 
radiation. 


1.  Introduction 

In  boundary  layer  phenomena,  radiation  and  chemical  effects  on  unsteady  MHD  fluid  flow,  heat 
and  mass  transfers  in  porous  media  is  significant  because  of  its  influence  exponentially  stretching 
porous  surface.  The  unsteady  MHD  fluid  flow  study  keeps  paving  ways  for  further  studies  due  to 
its  enormous  usage  in  every  facets  of  life.  The  investigation  such  as  this  can  be  modeled  and 
solved  experimentally,  analytically  or  numerically  depending  on  the  nature  of  the  modeled 
equations.  In  doing  this,  the  existing  models  are  modified  or  adjusted  to  reflect  the  current  trends 
that  are  germane  to  the  issues  being  discussed.  In  this  presentation,  efforts  are  made  to  reflect  the 
effects  of  some  unsteady  MHD  parameters  in  explaining  pollutions  in  porous  media.  MHD  fluid 
flow  in  porous  media  has  wide  and  numerous  applications  in  industry  and  environments  [1], 
Some  of  other  areas  of  applications  of  this  study  though  not  limited  are  the  flow  of  ground  water 
through  soil  and  rocks  (porous  media)  that  are  very  important  for  agriculture  and  pollution  control; 
extraction  of  oil  and  natural  gas  from  rocks  which  are  prominent  in  oil  and  gas  industries; 
functioning  of  tissues  in  body  (bone,  cartilage  and  muscle  and  so  on)  belong  to  porous  media,  flow 
of  blood  and  treatments  through  them;  understanding  various  medical  conditions  (such  as  tumor 
growth,  a  formation  of  porous  media)  and  their  treatment  (such  as  injection,  a  flow  through  porous 
media  in  medical  sciences). 

In  [1],  viscous  dissipative  heat  was  taken  into  account  under  the  influence  of  transverse  magnetic 
field,  [2]  reported  the  heat  and  mass  transfer  in  the  boundary  layers  on  an  exponentially  stretching 
continuous  surface.  Also,  [3]  presented  the  effects  of  heat  generation  and  thermal  radiation  on 
steady  MHD  flow  near  a  stagnation  point  on  a  stretching  sheet  in  porous  medium  and  presence  of 
variable  thermal  conductivity  and  mass  transfer,  it  was  discovered  that  temperature  increased  with 
increasing  radiation  parameter  R  and  concentration  decreased  with  increasing  Schmidt  number. 
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Unsteady  MHD  flow  and  heat  transfer  of  Nanofluid  over  a  permeable  shrinking  sheet  with  thermal 
radiation  and  chemical  reaction  was  analysed  [4].  The  analysis  of  unsteady  free  convective  flow 
and  mass  transfer  through  a  viscous,  incompressible  electrically  conducting  fluid  along  a  porous 
vertical  isothermal  non-conducting  uniformly  decaying  heat  generation  and  transverse  magnetic 
field  with  variable  suction  and  viscous  dissipation  was  performed  [5].  Unsteady  MHD  flow  of  a 
viscoelastic  fluid  along  vertical  porous  surface  with  fluctuating  temperature  and  concentration  was 
presented  [6].  [7]  verified  the  unsteady  MHD  flow  and  heat  transfer  for  Newtonian  fluids  over  an 
exponentially  stretching  sheet.  [8]  presented  thermal  radiation  effects  on  heat  and  mass  transfer  of 
MHD  flow  in  porous  media  over  exponentially-stretching  surface  but  neglected  unsteady  case.  In 
2018,  an  adaptive  MHD  convective  fluid  in  porous  media  was  presented  [9].  The  authors  critically 
discussed  the  applications  of  the  adaptive  model  to  environmental  issues.  In  view  of  the  above 
studies,  the  present  study  examined  the  radiation  and  chemical  effects  on  heat  and  mass  transfer  of 
unsteady  MHD  in  porous  media  over  an  exponentially  stretching  porous  surface.  Specifically,  the 
paper  contended  that  controls  of  environmental  pollution,  air  and  water  should  be  our  concern. 

2.  Materials  and  Method 

Consider  the  free  convective  thermal  radiation  effect  on  heat  and  mass  transfer  of  two  dimensional 
unsteady  MHD  flow  of  an  electrically  conducting,  unsteady,  viscous  and  incompressible  fluid  past 
an  exponentially  stretching  sheet  under  the  action  of  thermal  and  solutal  buoyancy  forces.  A 
uniform  transverse  variable  magnetic  field  B(x)  is  applied  perpendicular  to  the  direction  of  flow 
with  chemical  reaction  is  taking  place  in  the  fluid  flow.  The  flow  is  assumed  to  be  in  the  x- 
direction  with  y-axis  normal  to  it.  The  plate  is  maintained  at  the  temperature  and  species 
concentration  Tw,Cw  and  free  stream  temperature  and  species  concentration  TK ,  Cx  respectively. 

The  physical  model  and  equations  governing  the  radiation  and  chemical  effects  on  heat  and  mass 
transfer  of  two-dimensional  but  unsteady  MHD  fluid  flow  in  porous  media  are: 


Arrows  showing  B(x)  perpendicular  to  flow 


Figure  1  Showing  the  Physical  Model  and  coordinate  System 


—  +  —  =  0  ^ 

dx  dy 
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~r~  +  u  ~~  +  v  ~~  =  -(-oS(a-)  +  -^)w  +  v|-|-  +  gPT(T-Tx)+  gpc(c  -  Cj 
dt  dx  dy  p  K  dy" 


ar  ar  ar 

- h  U - h  V 

dt  8x  8y 


=  k^~8^-  +  Q„(T-Tj 

dy  dy 


dC  dC  dC  a2C  ,  X 

—  +  u  —  +  v —  =  D — -  -  y(C  -  C„  ) 
dt  dx  dy  dy" 

Subject  to  the  following  boundary  conditions: 


u  =  U0£  L,v  =  -V/L,T  =  Tw  =TX+  T0£ 2i ,  C  =  Cw  =  Cx  +  C0£ 2i  at  y  =  0 

DT  DC 

u  —>  0,  T  —>  T  ,  C  —>  C  as  y  qo,  —  ^  0,  —  ^  0,  —  ^  0 

dt  dt  dt 


(2) 

(3) 

(4) 

(5) 


where  u,  v,  C,  and  T  are  velocity  component  in  the  x  direction,  velocity  component  in  the  y 
direction,  concentration  of  the  fluid  species  and  fluid  temperature  respectively.  L  is  the  reference 
length  B(x)  is  the  magnetic  field  strength,  U0  is  the  reference  velocity,  V0  is  the  permeability  of 

the  porous  surface  respectively.  The  physical  quantities  K,  p,  v,  a,  D,  k,  C  ,  Q{)  and  y  are  the 

permeability  of  the  porous  medium  density,  fluid  kinematics,  viscosity,  electric  conductivity  of  the 
fluid,  coefficient  of  mass  diffusivity,  thermal  conductivity  of  the  fluid,  specific  heat,  rate  of 
specific  internal  heat  generation  or  absorption  and  reaction  rate  coefficient  respectively  while  g  is 

the  gravitational  acceleration,  /3T  and  fic  are  the  thermal  and  mass  expansion  coefficients 
respectively.  The  qr  is  the  radiative  heat  flux  in  the  y  direction.  By  using  the  Rosseland 
approximation  [3],  the  radiative  heat  flux  qr  is  given  by 

4cr0  dT 4  (6) 

Clr  = - 

3e>  dy 

Where  cr0  and  5  are  the  Stefan-Boltzmann  and  the  mean  absorption  coefficient  respectively.  We 

assume  the  temperature  differences  within  the  flow  are  sufficiently  small  such  that  T 4  may  be 
expressed  as  a  linear  function  of  temperature,  using  Taylor  series  to  expand  T 4  about  the  free 
stream  T x  and  neglecting  higher  order  terms,  this  gives  the  approximation 

T4=4T*T-3T*  (7) 


The  magnetic  field  B{x)  is  assumed  to  be  in  the  form  B(x)  =  B()£2L 

field  is  B() .  Introducing  the  stream  function  (//  (x,y)  such  that 

dip  dy/ 

u  = - ,  v  = - , 

dy  dx 


where  the  constant  magnetic 

(8) 


Substituting  (8)  in  (1),  continuity  equation  satisfies  Cauchy-Riemann  equations  and  equations 
(2)-(4),  gives 

dy/  ^dy/  d2y/  dy/ d2y/ _  a  ^  ^  dy/2  v  ^  dy/ 

dt  dy  dxdy  dx  dy2  p  °  i  J  k\  dy )  dy" 


v^Z  +  gPT(T-Tx)+gPc(C-Cx)  ^ 
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dy/  dy/  dT  dy/  dT 
dt  dy  dx  dx  dy 

a...  a...  a2. 


i6a0rj 


pc  P  ipcps  j 


dPT  ,  Q0 


dy2  pC , 


-(r-rj 


(10) 


dy/d2y/_  cr  ^ 

r dy/^ 

V 

2 dy/^ 

dx  dy2  p  ° 

K 

v  dy  p 

+  v  ^  +  gpT  {T  -  T ; )  +  g/?c  (c  -  c  J 
dy 


dyr  dy/  dT  dy/  dT 
dt  dy  dx  dx  dy 


i6a0rj 


pcP  ?>pcps 


~  P 
2 


dT+—(r-Tx) 


p~  j 


dy2  PC,, 


dC_  ^dy^d^  _dy^_dQ  =  Dd^  _  tc  c  \ 
dt  dy  dx  dx  dy  dy2  ^ 

The  corresponding  boundary  conditions  become: 

p  x  p  x  x  x 

-^  =  U0f\-^  =  V0f\T  =  Tw  =  TX+T0PL,C  =  Cw  =CX+  C0PL  aty  =  0 
dy  dx 

~~  ~pO,T  -a  Ty ,  C  -A  Cx  fli.yAO 
dy 


(11) 

(12) 


In  order  to  transform  the  equations  (9),  (10)  and  (11)  as  well  as  the  boundary  conditions  (12)  into 
an  ordinary  differential  equations,  [10]  introduced  the  following  similarity  transformations 
variables 


U  n 


y/{x,y)=pvU0U2Lf(r1),r1  =  y  UL£2L,T  =  Tx  +  To£^0(tj)C  =  C„  +  C0£ ^tj) 

V  2 vL 

The  equations  become 

+  2(f)2  -  (M  +  Da)f  +  Gr0  +  Gc(j)  =  0 


r 


\ 


1 +  -R 

V  3  , 


0'  +  Prf0'--Prf'0  +  PrQo0  =  0 


f  +  Seffi'--Sef't-SeW  =  0 

The  corresponding  boundary  conditions  take  the  form: 

/  =  /*,,  f  =  0-5,  0  =  0.5,  (j)  =  0.5,  atrj  =  0 
f  =  0,  0  =  0,  (j)  =  0  as  rj  — »  go 


(13) 

(14) 

(15) 

(16) 
(17) 


Where  u  is  unsteady  parameter,  where  M  is  the  Hartman  number,  Gr  is  the  thermal  Grashof 
number,  Gc  is  the  solutal  number,  Q  is  the  heat  source,  71  is  the  chemical  reaction  parameter,  Sc 
is  the  Schmidt  number,  Pr  is  Prandt  number,  Da  is  the  Darcy  porosity,  R  is  the  radiation 
parameter,  fw  is  the  permeability  at  the  plate.  The  equation  (14),  (15)  and  (16)  are  highly  non¬ 
linear  coupled  differential  equations  and  its  satisfying  the  boundary  conditions  (17)  are  solved 
numerically  by  applying  Nachtsheim-Swigert  shooting  iteration  technique  along  with  Runge- 
Kutta  fourth-order  integration  scheme. 
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3.  Results  and  Discussion 

Numerical  values  of  at  the  surface  for  different  values  of  Pr.  when  newly  introduced  parameter 
was  zero  is  presented  in  Table  1. 


Table  1:  Numerical  Comparison  of  the  present  study  at  the  surface  for  0’(0) 


Pr. 

M 

R 

-0'(o) 

Present  Study 

-0'(O)  [8] 

1 

0.08333 

0.09999 

2 

0.08333 

0.10000 

3 

0.08333 

0.10000 

1 

1 

0 

0.08333 

0.09999 

1 

0 

1 

0.08333 

0.09999 

1 

1 

0 

0.08333 

0.09999 

Therefore,  the  process  of  numerical  computation  for  the  skin  friction  coefficient,  the  Nusselt 
number  and  the  Sherwood  number  which  are  respectively  proportional  to  /"(()) ,  (7(0),  ^'(o) , 
from  this  study  at  the  plate  have  been  examined  for  different  values  of  the  parameters  are 
presented  in  Table  1  and  discussed. 


Table  2:  Effect  of  M,  G, ,  Gc,  Pr ,  Q ,  A,  Sc,  U,  R ,  Da ,  and  fw  on  /"( 0),  #'(o)  and  ^'(o) 


P 

VALUES 

F”(0) 

-0’(o) 

-0(0) 

M 

1 

-0.3893 

0.0728 

0.2165 

2 

-0.5647 

0.0376 

0.2096 

3 

-0.7167 

0.0089 

0.2042 

4 

-0.8515 

0.0150 

0.1998 

Gr 

1 

-0.2900 

0.0935 

0.2207 

2 

-0.0618 

0.1536 

0.2322 

3 

0.1474 

0.1927 

0.2411 

4 

0.3444 

0.2220 

0.2484 

Gc 

1 

-0.2900 

0.0935 

0.2207 

2 

-0.0781 

0.1349 

0.2296 

3 

0.1273 

0.1705 

0.2376 

4 

0.3271 

0.2016 

0.2457 

Pr 

0.71 

-0.2900 

0.0935 

0.2207 

3 

-0.3163 

0.2927 

0.2187 

3.5 

-0.3214 

0.3333 

0.2182 

3.9 

-0.3252 

0.3647 

0.2178 

Q 

-0.1 

-0.3297 

0.2582 

0.2145 

0 

-0.3198 

0.2143 

0.2160 

1 

-0.3002 

-0.0360 

0.2075 

1.1 

-0.2714 

-0.1578 

0.2108 

/l 

0.5 

-0.2900 

0.0935 

0.2207 

1 

-0.3006 

0.0857 

0.2828 

1.5 

-0.3086 

0.0804 

0.3336 

2 

-0.3152 

0.0765 

0.3774 

Sc 

0.24 

-0.2900 

0.0935 

0.2207 

0.62 

-0.3117 

0.0808 

0.3660 

0.78 

-0.3182 

0.0774 

0.4154 

2.62 

-0.3593 

0.0610 

0.8171 
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u 

i 

0.8293 

0.1544 

0.5941 

2 

0.9060 

0.0550 

0.5392 

3 

0.8911 

0.0470 

0.4792 

5 

0.8578 

0.0271 

0.3924 

R 

0.5 

0.8293 

0.1544 

0.5941 

1 

0.6443 

0.2283 

0.8088 

1.5 

0.3523 

0.2850 

1.0610 

2.5 

0.0734 

0.2777 

1.3388 

n 

1 

-0.2900 

0.0935 

0.2207 

3 

-0.0224 

0.2067 

0.2207 

5 

-0.9138 

0.0254 

0.1987 

7 

1.1373 

-0.0593 

0.9197 

f 

0.1 

-0.2900 

0.0935 

0.2207 

J  w 

0.5 

-0.3488 

0.1166 

0.2421 

1.0 

-0.4523 

0.1873 

0.2705 

1.5 

-0.6078 

0.3408 

0.3004 

Table  2  represents  the  numerical  results  of  variation  in  skin  friction,  Nusselt  and  Sherwood 
numbers  at  the  surface  with  M,  Da,  fw,  Sc,  Pr,  R  and  A  which  are  of  physical  and  engineering 
interest.  It  can  be  seen  from  the  results  that  an  increase  in  the  values  of  M,  Sc,  Pr,  Da,  fw,  decrease 
the  flow  boundary  layer  while  increase  in  values  of  Gr,  Gc,  Q,  A  and  R  increase  the  flow 
boundary  layer.  Table  1  depicts  that  an  increase  in  the  values  of  M,  Q,  Sc,  Da,  R  thicken  the 
thermal  boundary  layer  by  reducing  the  rate  at  which  heat  diffuse  out  of  the  system  while  increase 
in  Gr,  Gc,  Pr,  fw  and  Pr  reduces  the  thickness  of  the  thermal  boundary  layer.  Also,  the  results 
showed  that  increase  in  Gr,  Gc,  Q,  Sc,  fw,  R  causes  thinning  in  the  concentration  boundary  layer 
while  M,  Pr,  Da  thicken  the  mass  boundary  layer.  The  following  graphs  (Figures  2  -7)  further 
elucidate  the  research  work. 

O 

O  .-4- 

0.3 

O  _2 


O  .  1 

O 

0  12  3  4  5  6 


Figure  2:  Velocity  profile  for  the  variation  of  Hatmann  number  M 
Figure  2  shows  that  an  increase  in  Hartmann’s  number  led  to  decrease  in  velocity  profile 


Figure  3:  Velocity  profile  for  the  variation  of  Darcy  Porosity  Da 


Figure  3  shows  that  an  increase  in  Darcy  porosity  number  led  to  decrease  in  velocity  profile 
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Figure  4:  Temperature  profile  for  the  variation  of  the  Chemical  reaction  2 

Figure  4  shows  that  an  increase  in  chemical  reaction  Tied  to  increase  in  temperature  profile.  Its 
influence  thereby  increasing  the  rate  at  which  the  fluid  was  able  to  flow  out  of  the  system. 


Figure  5:  Temperature  profile  for  the  variation  of  the  Prandtl  number  Pr. 

Figure  5  shows  that  an  increase  in  Prandtl  number  Pr.  led  to  decrease  in  temperature  profile.  In 
this  case  Pr.  decrease  the  thickness  of  thermal  boundary  layer  and  as  a  result,  heat  was  able  to 
diffuse  out  of  the  system  thereby  inceasing  the  temperatue  profile. 


Figure  6:  Temperature  profile  for  the  variation  of  the  Darcy  Porosity  Da 
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Figure  6  shows  that  an  increase  in  Darcy  Porosity  Da.  led  to  increase  in  temperature  profile.  Its 
effect  on  hydrological  flow  in  sand  soil  did  not  permit  turbulence  hence  porosity  influences 
laminar  flow. 


Figure  7 :  Concentration  profile  for  the  variation  of  the  Radiation  R 

Figure  7  shows  that  an  increase  in  Radiation  R.  led  to  decrease  in  concentration  profile.  Its 
decrease  causedthe  thinning  of  mass  boundary  layer 

4.  Conclusion 

In  this  study,  a  numerical  analysis  was  presented  to  investigate  the  influence  of  radiation  and 
chemical  reaction  effects  on  unsteady  MHD  boundary  layer  fluid  flow  in  porous  media.  The  non¬ 
linear  and  coupled  governing  equations  were  solved  numerically  by  using  fourth  order  Runge- 
Kutta  scheme  using  Maple  18  software.  Velocity,  temperature  and  concentration  profiles  were 
presented  graphically  and  explained  appropriately.  The  fundamental  independent  parameters 
found  to  effect  the  problem  under  consideration  were  the  Magnetic  field  parameter,  Darcy 
porosity  parameter,  chemical  reaction  and  Radiation  parameter. 
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